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Abstract 

Events with tagged photons in the process of electron-positron annihilation into 

hadrons are considered. The initial state radiation is suggested to scan the hadronic 

cross section with the energy. QED radiative corrections are taken into account. The 

■^j- ■ results for the total and exclusive cross sections are given in an analytic form. Some 

■ numerical estimates are presented. 
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> '. 1 Introduction 

Experiments with tagged photons, radiated from the initial state in electron-proton and 
electron-positron collisions, can become particularly attractive. The reason is that these 
radiative processes will permit to extract information about the final states at continuously 
varying values of the collision energy. To investigate deep inelastic scattering the authors of 
Ref. U suggested to use radiative events instead of running colliders at reduced beam energies. 
The method takes advantage of a photon detector (PD) placed in the very forward direction, 
as seen from the incoming electron beam. The effective beam energy, for each radiative event, 
is determined by the energy of the hard photon is observed in PD. In fact, radiative events 
were already used to measure the structure function F 2 down to Q 2 > 1.5 GeV 2 [0]. The 
specific theoretical work concerns the evaluation of QED radiative corrections |3], Q to the 
radiative Born cross section. With an accurate determination of the cross sections and of 
the possible sources of background we believe that the use of radiative events may become 
particularly useful to carry investigations at various present and future machines. 

The important role of the initial state radiation in the process of electron-positron an- 
nihilation was underlined in a series of papers by V.N. Baier and VA. Khoze 0], where 
the radiative process was studied in detail in the Born approximation. In these papers the 
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mechanism of returning to a resonant region was discovered. This mechanism consists in the 
preferable emission of photons from the initial particles, which provides a resonant kinematics 
of a subprocess. A utilization of radiative events can become a common type of investigations 
at various machines. 

In this paper we derive explicit formulae for the spectrum of tagged photons. The calculati- 
ons are performed having an accuracy of the per-mille order as an aim. Formulae can be used 
at electron-positron colliders to investigate, for instance, hadronic final states at intermediate 
energies. A measurement of the total hadronic cross section at low energies is essential for 
an high precision test of the Standard Model particularly for a precise determination of the 
fine structure constant cxqed^Mz) and of the muon anomalous magnetic moment (g — 2) M . 
The largest contribution to the errors for these quantities comes from the large indetermina- 
tion still present on the measurement of the total hadronic cross section in electron-positron 
annihilation at center-of-mass energies of a few GeVs. We will consider here the radiatively 
corrected cross section for the electron-positron annihilation process 

e _ (pi) + e + (p 2 ) — > 7 (A;) + H(q), k = (1 - z) Pl , (1) 

where if is a generic hadronic state. The hard photon hitting the photon detector has a 
momentum k and an energy fraction 1 — z with respect to the beam energy. In the following 
we assume that the photon detector is placed along the electron beam direction, and has an 
opening angle 26$ <C 1, such that e 2 9 2 ^> m 2 , with m the electron mass, and e the beam energy. 
To evaluate the process with an accuracy of the per mille requires a careful investigation of the 
radiative corrections. This paper is organized as follows. In Section 2 we consider the cross 
section of the process ([!]) in the Born approximation. We give formulae suitable to study as 
differential distributions in hadronic channels, as well as the total (in terms of quantity R) and 
inclusive (in terms of hadron fragmentation functions) hadronic cross sections. In Sec. 3 we 
calculate separate contributions into radiatively corrected cross section of process (H) within 
the next-to-leading accuracy. In Sec. 3.1 the contribution due to virtual and soft photon 
emission is investigated. In Sec. 3.2 the case, when additional hard photon hits a photon 
detector is considered. In Sec. 3.3 the contribution due to hard photon emission, which does 
not hit a photon detector is derived. In Sec. 4 we sum up all the contributions and give the 
final result. In Conclusions we summarize the results and give some numerical illustrations. 

2 The Born approximation 

In order to obtain the Born approximation for the cross section of the process ([!]), when the 
PD is placed in front of electron (positron) beam, we can use the quasireal electron method |6j . 
It gives 

da(k,p 1 ,p 2 ) = dW pi (k)a (p 1 - k,p 2 ), (2) 

where dW Pl (k) is the probability to radiate photon with energy fraction 1 — z inside a narrow 
cone with the polar angle not exceeding fl <C 1 around the incoming electron, and dcr is 
the differential cross section for the radiationless process of electron-positron annihilation into 
hadrons at the reduced electron beam energy. The form of both, dW Pl (k) and (?o(pi — k,p 2 ) 
is well known: 

a 1 + z 2 2z e 2 9 2 

dW pl (k) = —P 1 (z,L )dz } Pi ( Z)Lo ) = -^ Lq - , Lo = i n ^£. ( 3 ) 
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We need further the general form of the lowest order cross section o" for the process e + (ziP2) + 
e~(zpi) —> hadrons boosted along the beam axis (pj: 

2 ~ W<y ~ PiQ 

q = zp± + z x p 2 q = sz x z, g pa = g pa - , p lp = p lp - —q p , 

where q is the full 4-momentum of final hadrons, qj is 4-momentum of an individual hadron, 
s = 2pip 2 = 4:E 2 is the full center-of-mass energy squared, and H pa is the hadronic tensor. 
The vacuum polarization operator II(g 2 ) of the virtual photon with momentum q is a known 
function [0] and will not be specified here. 

The tensors H pa and L pa obey the current conservation conditions once saturated with the 
4-vector q. The differential cross section with respect to the tagged photon energy fraction z 
can be obtained by performing the integration on the hadrons phase space. It takes the form 

^=^ Pl ( z ,L ) * (z,l). (6) 

0.Z Z7T 

Each hadronic state is described by its own hadronic tensor. The cross section in Eqs. (3) 
and (4) is suitable for different uses and, as mentioned above, it can be used to check different 
theoretical predictions. 

The sum of the contributions of all hadronic channels by means of the relation 

[H pa dr = f h (q 2 )g P «, (7) 

h J 

can be expressed in terms of the ratio of the total cross section for annihilation into hadrons 
and muons R = Ohjo^. For the final state we get 



and so, 

Afe 2 ) = ^f 1 ^ 2 )- (8) 

Substituting this expression into the right hand side of Eqs. (3,4) results in the replacement 
a (z,z 1 )=R(q 2 )A7ra 2 K(q 2 )/(3q 2 ). 

In experiments of semi-inclusive type one fixes an hadron with 3-momentum q 1 energy ei 
and mass M in every event and sum over all the rest. In this case instead of Eq. (|7|) we will 
have (similarly to the Deep Inelastic Scattering (DIS) case J|]): 



V f H dT = d3gi (9) 

ZfJ pa pa 2 £i (2tt) 3 ' [ ) 

^F 2 ( V ,q)q lp q lp , r, = — 



4 q 2 

HpJ = F i(v, q 2 )g P * - — F 2(v> q 2 )qi p qi p , v = t^t > 1 
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where we have introduced two dimentionless functions Fi(r), q 2 ) and -F^??, Q 2 ) in a way similar 
to the DIS case. 

By introducing the dimentionless variable A = 2gg 1 /(2zp 1 g 1 ), we can write the correspon- 
ding cross section for radiative events in e + e~ annihilation in the same form as in the case of 
deep inelastic scattering with a tagged photon ||: 

do" a 2 (q 2 ) a 2 1 d 3 q 1 

Tz = ^r^ Pi(z ' Lo)s( ^ A ' g) w^r' (10) 

r( \ *\ j-i / 2 , , 2F 2 {r 1 ,q 2 ) ( x M 2 2 2 \ 

Z{r),\,q) = F 1 ( n ,q)+ [\ - 1 - —r] A j. 



3 Radiative corrections 

For the radiative corrections (RC) to the cross section (5) we will restrict ourselves only to 
terms containing second and first powers of large logarithms L, and omit terms which don't 
contain them i.e. we will keep leading and next-to-leading logarithmic contributions. We will 
consider in subsection 3.1 the contribution from one-loop virtual photon as well as from the 
emission of soft real ones. In 3.2 we will discuss the double hard photon emission process. 



3.1 Corrections due to the virtual and real soft photons 

The interference of Born and one-loop contributions to the amplitude of the initial state radia- 
tion in annihilation of e + e~ into hadrons can be obtained from the analogous quantity of hard 
photon emission in electron-proton scattering [|J. We do that by using the crossing transfor- 
mation. For the contribution coming from the emission of real soft photons a straightforward 
calculation gives: 



da s a 



dan n 



mAe 1 9 7r 2 

2(L S - I n— — + -L 2 

v s ' Xe 2 s 3 



s 4 
In — = L + L g , L e = ln — . (11) 



where A is the photon mass, Ae is the energy carried by the soft photon. The sum of the two 
contributions is free from infrared singularities. It reads 



where 



7T 



2 



9 Ae 



p = 4(L s -l)lnA + 3L g - — --, L q = L s + \nz, A = — < 1, (13) 

where k and u are the 3-momentum and the energy of the hard photon respectively. The ten- 
sors A pa and B pa have a rather involved form. The first can be obtained from the corresponding 
expressions of Ref. ||. The tensor B pa coincides with the one of the Born approximation. In 
the kinematical region where the hard photon is emitted close to the initial electron direction 
of motion one has 



\ _|_ z 2 2m 2 z\ 2 
Bpa = ~z \ yi (l -z) ~ ~yf) L ^ q ^ A p° = ^ A 9 L pM)i q = z Pl +p 2 , (14) 
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where tensor L pa is given in Eq. (4), yi = 2kp±, and quantity A g reads 



Azsm 2 s 
5 — L s lnz H 

Vi Vi 



1 + z 2 



+ 2 in — laz ) + 



2 1 



;-2L s In z - In 2 z + 2Li 2 (l - z) 

Li 2W = -/'d t "' (1 - tx) . 

Jo t 



(15) 



Further integration over the hard photon phase space can be performed within the logarithmic 
accuracy by using the integrals 



d 3 k 
2ixk n 



1 m 2 ln(yi/m 2 ) 



' 2 ' 

yi m 



L l i r2 i r W1 

— Ln, — ; r, — L n H Ln In 1 

2 °' 2(1 -z '4 2 v 



dz. 
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The final expression for the Born cross section corrected for the emission of soft and virtual 
photons has the form 



da B+V+S 



dz 



N 



a (z, 1) 
1 + z 2 



a 
2^ 



Pi(z,L c 



a V 

IT 1 



N) 



+ 



1 - z 
1 + 2z - z 2 

2(1 -z) 



(L + ln z) ln z - — + 2U 2 (z) 



L -2P 1 (z,L )\n^ 



Az 

L + L Q \nz. 

1 — z 



(17) 



3.2 Two hard photons are tagged by the detector 

If an additional hard photon emitted by the initial-state electron hits the PD, we cannot use 
the quasireal electron method and have to calculate the corresponding contribution starting 
from Feynman diagrams. 

We can use double hard photon spectra as given in Ref. |J for annihilation diagrams only 
and write the cross section under consideration as follows 



do| 

dz 



*o(*,l)(£) L 



l-z-A 



dx 

T 



7T 



[z A + (1 



X 



In 



(1 -x) 2 (l -z-x) 



zx 



+ zx(l 



x) — x 2 (l — x — z) A 



t 



x(l — x) (1 — z — x) 
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2r(l 



:i-x) 



x 



;is) 



z 2 + a-x) 2 . 



Here the variable x under the integral sign is the energy fraction of one hard photon. The 
quantity 1 — z — x is the energy fraction of the second hard photon provided that their total 
energy fraction equals to 1 — z. We write the index cl in the left hand side of Eq. flT8| ) to 
emphasize that this contribution arises from the collinear kinematics, when the additional 
hard photon is emitted along the initial electron with 4-momentum p±. 
The integration in the right hand side of Eq. (pi) leads to the result 



dof 

dz 



a (z, 1) 



a 



■L 



2nJ 2 



PS\z)+2 



1 + z 2 



1 - z 



+ 6(1 -z) + 



3 + z 2 
1 - z 



\n 2 z 



4(1 + z) 



In 



ln z 



1 - z 



21nA 



A 



(19) 



(2) 

where the quantity P@ (z) represents the so called O-term of the second-order electron struc- 
ture function: 



1 + z 2 
' 1-z 



In 



— + + lnz-2(l-z). 

z 2) 



(20) 



3.3 Additional hard photon is emitted outside PD 



If an additional hard photon, emitted from the initial state, does not hit the PD situated in 
the direction of motion of the initial electron we distinguish the case when it is emitted in the 
direction close, within a small cone with angle 9' <C 1, to the direction of the initial positron. 
In this case we obtain: 



dz 



ry pl — S/z ry 

= —P^z, Lq) J a — Px(l - x, L')a (z, 1 - x)dx, 



(21) 



where L' = L s + ln(# /2 /4), 5 = M 2 /s, and M 2 is the minimal hadron mass squared. We 
suppose that z ~ 1. 

We have introduced the additional auxiliary parameter 0' < 1 which, together with 9 , 
separates collinear and semi-collinear kinematics of the second hard photon. Contrary to 9q, 
which is supposed to determine the PD acceptance, 9' will disappear in the sum of the collinear 
and semi-collinear contributions of the second photon. This last kinematical region gives 



da 



H 



a ) 2 P(z 1 ) f 167r2ft2 T(r z r) 

H pa (q 2 )dT(q 2 ) 



dz 

T(c,z,x) 



1 



-((z-x 2 ) 2 + z 2 (l-x 1 ) 2 )g prT 



s(q 2 2 ) 2 \l-U(q 2 )\ 2 
+ 2(z(l - x x ) - x 2 ) (z 2 p lp p la + p 2p p 2a ) 
q 2 = zp 1 +p 2 -k 1 , c = cosk 1 p 1 . 
The phase volume of the second photon is parametrized as: 



x 



x 



zi = o(l-c), x 2 = -(l + c), 



(22) 



dx 



2-k 



1-61/2 



27TU0 3 J X J 2ix 

A _i + 0/2/ 2 



dc, 



x = 



2(z-S) 



l + z + c(l-z)' 



(23) 



Explicitly extracting the angular singularities we represent this expression as 

, z-6 



da 



H 



dz 



T, sc (z) + In 



91 I 



dx z 2 + (z — x\ 



, . x 



z z 



a (z - x, 1) 



l-S/z 

+ ln V 2 J ^a + a-^VoM-*) 

A 



(24) 



87r 2 a 2 



dc 



-i 



drr 



x 



T(c, z, x) - T(l, z, x) T(c, z, x) - T(-l, z, x) 



+ 



1 + c 
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4 Complete QED correction and leading logarithmic 
approximation 

The final result in the order 0(a) for radiative corrections to radiative events can be written 
as follows: 



dcr a, „ . . . . . 

— = —P x {z,L Q )a Q {z,l){l + r) 



(25) 



^-P x {z, L )a (z, 1) + ^ 2 ^L (h pW(z) + Gy o (z, 1) 

^j" Ci(x)a a {z,l - x)Ax + L„ 'J C 2 {z, x)cr„(z - x, l)dx + £, 



+ Pi(z,L. 



where the last term is defined by Eq. (23) and 



C x {x) = P 1 (l-x,L s )Q(x-A) + (L s -l)(2\nA + -)8(x), 



C 2 (z,x) 
G(z) 



z 2 + (z — x) 2 

Z 2 X 



G(x - A) + (2 In A + - - 2ln z)8(x), 

= l±i!( 3 \nz - 2Li 2 (z)) + -(1 + z) In 2 z - 2(1 + z) 
l-z K y T ' 1-z 

1 — 16z — z 2 4z\nz 
+ —T7z + 



(26) 



ln(l - z) 



2(1 -z) 



1 - z 



In order to include the higher order leading corrections to the tagged photon differen- 
tial cross-section and show the agreement of our calculation with the well known Drell-Yan 
representation for the total hadronic cross-section at electron-positron annihilation [110] 



i i 



18} 



J dxi J dx 2 D(xi, a c s)D(x 2 , a e s)a(xix 2 s) 



(27) 



5 S/xi 

where the electron structure functions include both nonsinglet and singlet parts 

D(xt, a cS ) = D NS (x, a cS ) + D s (xi, a cS ), 

it is convenient to introduce the quantity 

i i 

S = D(z, a cS ) J dxi J dx 2 D(x 1 ,a cS )D(x 2 ,a cS )a (zx 1 ,x 2 ). 

S/z 5 j zx i 

Note that the shifted cross-section in Eq. ([27]) has just the same meaning as in Eq. 
a(x l x 2 s) = a (xx,x 2 ). 

The structure functions [11], |12| , which enter into the right side of Eq. ( p8[) , are 



(28) 



(29) 



n=l 



(30) 
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2 1 / \ 3 



2Pi ® P(x) - -R(x) 



(31) 



where 



P x (x) = \im\-±^e(l-A-x)+(~+2\nA)8(l-x)\, 
a^o L 1 — x V2 / J 

P(x) = 2(l + x)lnx + ^— -(4 + 7x + 4x 2 ), 

i 

Pf" = P 1 (x)^---0P 1 (x) , Pi(x)®P l (x)=Jp 1 (t)P 1 (j)j 



and the effective electromagnetic couplings in the right side of Eq. ( ^9|) are 

tteflf = -37T 111 (l - 7^0^ , «eff = -37T In Q - j^" ) , «eflf = -37T In ^1 - • (32) 

At fixed values of z (z < 1) the quantity E defines the leading logarithmic contributions 
into differential cross-section for the events with tagged particles. That corresponds to only 0- 
terms in the expansion of the structure function D(z, a e g) before the integral sign in Eq. (|29|) . 
If we consider photonic corrections (as in previous Sections) it needs to restrict ourselves with 
the nonsinglet part of the electron structure functions and with the first order terms in the 
expansion of all effective couplings, namely: 

a eS — > aL , a cS — > aL e , a eS — > aL s . (33) 



It is easy to see that in this case the leading contribution into differential cross-section (|25| ) 
can be obtained as an expansion of the quantity £(z < 1) by the powers of a, keeping the 
terms of the order a 2 in the production of D-functions. 

If we want to include the contribution due to e + e~-pair (real and virtual) production it 
is required to use both nonsinglet and singlet structure functions and effective couplings 
defined by Eq. ([32]). Note that the insertion into consideration of higher order corrections rises 
additional questions about concrete experimental conditions concerning registration of events 
with e + e~-pairs. 

The total hadronic cross-section in e + e~ annihilation can be obtained by integration of 
quantity £ over z 

i 11 

J dz D(z,a c s) J dxi J dx2 D(x\, a c g)D(x2, a c s)a(zxiX2s). (34) 

S 5/z 8/zxi 



cr(s 



We can integrate the expression in the right side of Eq. (p4|) over the variable z provided the 
quantity zx\ = y is fixed 

i i 11 

J dz D(z, a cS ) J dxt D(x u a cS ) = J dz J dy D(z, a eS )D 5 eff j (35) 

S 8/z 8 y 

1 

= J dy D(y, a cS + 5 eff ), a eS + a eS = a eS . 
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Using this result and definition of a e g we indicate the equivalence of the Drell-Yan form of 
the total cross-section as given by Eq. (|27| ) and the representation of the cross-section by 
Eq. (0). 

Let us show now that D-functions in expression for the quantity £ have effective couplings 
as given by Eq. (|3"2|). By definition the nonsinglet electron structure function satisfies the 
equation [K| 

s 1 
1 f dsi r dz Si 



D(x,s,s )=5(l-x) + — / — a( Sl ) / — D(z)D( — , — ), (36) 

ZlX J S\ J Z VZ So- 



so x 



where a(si) is the electromagnetic running coupling 

a(si) =a(l - ^-ln^r 
V Sir m z 



-i 



and so(s) is the minimal (maximal) virtuality of the particle, which radiate photons and 
e + e~-pairs. 

The structure function D(z, a e g) describes the photon emission and pair production inside 
narrow cone along the electron beam direction. In this kinematics Sq = m 2 , s = e 2 8q. The 
corresponding iterative solution of the Eq. fl3"6|) has the form (|30| ) with a e g = a e g. The 
structure function D(xi,a e g) describes the events, when emitted (by the electron) particles 
escape this narrow cone. In this case Sq = e 2 8 2 ., s = 4e 2 . The corresponding solution 



of Eq. (|36"D gives the structure function with a e s = a e g. At last, the structure function 
D(x2, a e ff) is responsible for the radiation off the positron into the whole phase space. In this 
case s = m 2 , s = As 2 . Therefore we obtain D-function with a e s = a e fr- The analogous 
consideration can be performed for the singlet part of structure functions. 

When writing the representation (0) for the total cross-section we, in fact, divide the phase 
space of the particles emitted by the electron on the regions inside and outside the narrow 
cone along electron beam direction. Therefore we can use this representation to investigate the 
events with tagged particles in both this regions. As we saw before the differential cross-section 
for events with tagged particles inside the narrow cone is defined by the quantity E(z < 1). 
In order to obtain the corresponding differential cross-section for events with tagged particles 
outside this narrow cone we have to change the places of a e g and 5 e ff i n expression for £(z, 1). 
This follows from the symmetry of representation (Rl) relative such change. 



5 Conclusion 



In sum the formulae (|35| , |29"D are the main results of our paper. In Fig. 1 we show the cross 
section da/dz as a function of z. The beam energy is chosen to be -Ebeam = 0.5 GeV. The 
region of z values is limited by the pion production threshold at the left, and by the threshold 
of photon detection (we choose 50 MeV) at the right. The peak in middle corresponds to the 
large contribution of the p-meson. Values of R, used for numerical estimations, were taken 
from [15]. The values of corrections r (see Eq. (|35|) ) in percent are shown in Fig. 2. 



So, we calculated the cross-section of e + e annihilation with detection of a hard photon 
at small angles in respect to the electron beam. The general structure of a measured cross- 
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Figure 1: The cross section of e + e — ► hadrons with tagged photon, 
section, from which one should extract the annihilation crossection ctq, looks as follows: 



a 



a 1 -L + b 1 - + a 2 - L + & 2 1 — I £ + 



7T 



7T 



7T 



7T 



7T 



(37) 



where L denotes some large logarithm. We calculated the terms ax, b\, a 2 , &2 and some 
contributions to c 2 . The generalized formula ( p9l) allows to involve the leading terms of the 
order 0(a 3 L 3 ). In this way our formulae provide high theoretical precision. 

Similar formulae can be obtained for an experimental set-up by tagging a definite hadron. 
By using e+e--machines such as BEPS, DA$NE pj, VEPP, CLEO, SLAC-B/factory and 



others with luminosities of order 10 cm s~ , one may be in principle able to scan, by mea- 
suring the initial state radiation spectrum, the whole energy region of hadron production with 
an effective luminosity of the order of 10 31 cm 2 s _1 . We hope further study will follow on these 
issues both from the experiments and theory. 
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